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We have created one-, two-, and three-dimensional quantum gases and study 
the superfluid to Mott insulator transition. Measurements of the transi- 
tion using Bragg spectroscopy show that the excitation spectra of the low- 
dimensional superfluids differ significantly from the three-dimensional case. 
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1. LOW-DIMENSIONAL QUANTUM SYSTEMS 

A low dimensional gas can be created in a trap when the confining po- 
tential restricts the motion of the particles to one or two dimensions with 
the other motional degrees of freedom being frozen out. Quantum mechan- 
ically this implies that the particles occupy only the ground state in the 
restricted directions, i.e. both, the thermal energy ksT and the interac- 
tion energy \x have to be much smaller than the energy level spacing. In 
general, low-dimensional systems exhibit increased quantum fluctuations of 
the phase. For a homogeneous 2D gas these phase fluctuations allow Bose- 
Einstein condensation only at zero temperature and for a homogeneous ID 
gas no Bose condensation is possible at all. For trapped low-dimensional 
gases the Bose-Einstein phase transition takes place at finite temperatures 
for both, and but the fluctuating phase alters the properties of 
the gas. Low-dimensional quantum systems exhibit a wealth of fascinating 
phenomena whose explanations go beyond the mean-field descriptiorP. For 
example, in one dimension bosons acquire fermionic properties at low den- 
sities and strong interactions, forming the so-called Tonks- Girardeau gas 4 . 
For two dimensions a Berezinskii-Kosterlitz-Thouless transition is predicted 
resulting in the spontaneous formation of vortices from thermal excitations^. 
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A one-dimensional gas can be realized in a cigar shaped harmonic trap- 
ping geometry characterized by the frequencies the two strongly con- 
fining axes and the weakly confining axis. A two-dimensional trap 
can be realized in a disk shaped geometry with one strongly confining axis 
characterized by ui 2D and two weakly confined directions by uPf* . The con- 
ditions for achieving the low-dimensional quantum regime then read 



, „ r fiuj] for ID; 

kBT >^\hJ D for 2D. (1) 



A one-dimensional trapped Bose gas in the weakly interacting regime was 
recently created and studiecP. This experiment revealed the distinctively 
different excitation spectrum of a one-dimensional quantum system as com- 
pared to three dimensions. Previous studies accessed a regime, where a 
Bose condensate with \x < hco]® coexisted with a three-dimensional thermal 
cloud^. Also Bose gases were created in optical lattices with significant tun- 
nelling between the one-dimensional subsystems^. Theoretical understand- 
ing of the homogeneous one-dimensional Bose gas was pioneered by Lieb and 
Liniger in the 1960s^ and recently adapted to atom trapa^ ^ 2 * 13 * 14 * 15 ^. A 
fundamental feature of the one-dimensional Bose gas is that its ground state 
energy and excitation spectrum are determined by a single parameter 7 over 
the whole range from weak to strong interactions. 7 = m 9lD is the ratio 
between interaction energy and kinetic energy (m: atomic mass, g\D'- ID 
coupling constant, n\D\ fD density). The transition from weak to strong in- 
teractions can be studied by applying a lattice potential along the symmetry 
axis of the ID gasP. 

Quasi-homogeneous two-dimensional quantum gases have been studied 
for several years using thin films of liquid Helium on a surface^. With 
ultracold alkali atoms quantum degenerate atomic gases in two dimensions 
have been realized in a standing wave laser fielcPl Moreover, the regime with 
< fiu!^ D has been studied using 23 Na in a disk shaped optical trap where 
features of two-dimensional condensate expansion were observed^. With 
133 Cs a quantum degenerate sample has been prepared in an optical trap 
close to a surface^ and in extremely fast rotating condensates the two- 
dimensional regime has also been approache <P. In this paper we present 
experiments with Bose gases fully in the 2D regime with ksT '/Huj 2D < 6 x 
1(T 3 and fi/huj 2 z D < 0.1. The gas is subject to a periodic potential along its 
weakly confined directions, which allows us to tune the gas from a weakly 
interacting superfluid to a strongly correlated Mott insulator. 
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2. SUPERFLUID TO MOTT INSULATOR TRANSITION 

Bosons which are localized in the minima of a periodic potential may 
be described by the Bose-Hubbard Hamiltoniad 22 * 23 ! 

H = -J^ajdj + Y e i n i + ~2 Y ni ( ni ~ 1 )- ( 2 ) 

i,j i i 

J denotes the hopping matrix element between neighboring lattice sites that 
determines the rate of which a particle disappears from lattice site i and 
tunnels to the adjacent lattice site j (Si and a\ are the annihilation and 
creation operators for an atom at lattice site i, respectively). The total 
tunnelling rate J includes the possibly anisotropic tunnelling to all next 
neighboring lattice sites: J = 2{J X + J y + J z ). ej characterizes the inhomo- 
geneity of the atom trap, is the occupation number of lattice site i and 
U is the onsite-interaction energy between two bosons at the same lattice 
site. The Bose-Hubbard Hamiltonian exhibits a quantum phase transition 
between the superfluid state and a Mott insulator for a critical ratio (U/J) c 
which has been observed in an optical lattice, where the periodic potential 
is created by standing wave laser 

fieldiP. 

The intensity of the laser fields 
determines the value of U/J. According to a mean-field calculation one 
expects (U/J) c = 5. £p 2 l 23 ! This approximation is valid for two and three 
dimensions, but for one spatial dimension stronger quantum fluctuations 
of the phase must be taken into accounlP^. Quantum Monte-Carlo simu- 
lations showS that these lower the onset of the Mott insulating phase to 
(U/J) c — 1.8. We have recently obtained experimental evidence for these 
fluctuations^. 
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Fig. 1. Quantum phase transition from a superfluid to a Mott insulator in 
three dimensions. The absorption images are taken after 25 ms of ballistic 
expansion. The interference peaks in a) are separated by a momentum of 
2hk, where k = 2ir/\ is the wave vector of the laser creating the optical 
lattice. 
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One signature of the phase transition can be observed by studying bal- 
listic expansion of the atoms from the lattice. For strong tunnelling, the 
system is superfluid, and the wave function of an atom is delocalized over 
many lattice sites. The absorption image shows interference peaks with 
a periodicity according to the lattice vector (see Figure^,). In contrast, 
for strong interactions in the Mott insulating phase the atoms are localized 
within a single potential well and have no definite phase relation to the 
next potential well. Therefore, no interference pattern develops during the 
expansion (see Figure^). 



3. EXPERIMENTAL SETUP 
3.1. Bose condensates in an optical lattice 

In the experiment, we collect up to 2 x 10 9 87 Rb atoms in a vapor cell 
magneto-optical trap. After polarization gradient cooling and optical pump- 
ing into the \F = 2,mf = 2) hyperfine ground state the atoms are captured 
in a magnetic quadrupole trap. After magnetic transport of the trapped 
atoms over a distance of 40 cm the magnetic trapping potential is converted 
into the harmonic and elongated potential of a QUIC trapP^. Subsequently, 
we perform radio frequency induced evaporation of the cloud over a period 
of 25 s. After evaporation we observe almost pure Bose-Einstein condensates 
of up to 1.5 x 10 5 atoms. Following the condensation we adiabatically change 
the trapping geometry to an approximately spherical symmetry with trap- 
ping frequencies of u x = 2ir X 18 Hz, uj y = 2ir X 20 Hz, and lo z = 2ir x 22 Hz. 
This reduces the peak density by a factor 4 and allows us to load the optical 
lattice more uniformly. 

The optical lattice is formed by three retro-reflected laser beams. Each 
beam is derived from a laser diode at a wavelength of A = 826 nm. At the 
position of the condensate the beams are circularly focused to l/e 2 -radii of 
120 //m (x and y axes) and 105 /iin (z). The three beams possess mutually 
orthogonal polarizations and their frequencies are offset with respect to each 
other by several ten MHz. We stabilize the lasers to a high- finesse Fabry- 
Perot cavity, thereby reducing their line width to ~ 10 kHz. In order to load 
the condensate into the ground state of the optical lattice, the intensities 
of the lasers are slowly increased to their final values using an exponential 
ramp with a time constant of 25 ms and a duration of 100 ms (see Fig. E}. 
The resulting optical potential depths V x ,y,z are proportional to the laser 
intensities and are conveniently expressed in terms of the recoil energy Er = 

with k = ^r- and the atomic mass m. 
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3.2. Preparation of low-dimensional quantum gases 

Using the optical lattice we realize one-, two-, and three-dimensional 
quantum gases. To prepare an array of one-dimensional tubes, two lattice 
axes are ramped to a fixed value V± = V x = V z = 30 En. In this config- 
uration the transverse tunnelling rates J x and J z are small and contribute 
only a correction of the order of J x ,z/ M C 1 to the ID characteristics of 
the individual tubes, where [i is the chemical potential of the sample. To 
study the superfluid to Mott insulator transition in this system we apply 
an additional optical lattice along the symmetry axis with a potential depth 
V y <C V±. Upon variation of V y the state of the system changes from super- 
fluid to Mott insulating. Experimentally we observe the same behavior of the 
system for the following preparation schemes: (i) when ramping all lattice 
laser beams at the same time, (ii) when first preparing a one-dimensional gas 
and subsequently ramping up V y and (iii) when creating a three-dimensional 
lattice first and subsequently lower V y to the desired value. 

We prepare two-dimensional quantum gases by applying one strong lat- 
tice laser (V z = 30 En) to freeze the atomic motion in the z-direction only. 
The two other laser beams are kept at lower levels with V x = V y . For the 
three-dimensional situation we use the same power in all three laser beams 

(V x =Vy = V Z ). 

4. BRAGG SPECTROSCOPY IN AN OPTICAL LATTICE 
4.1. Principle of the measurement 

The excitation spectrum reveals fundamental properties of the quan- 
tum state of the s ystem . We study the excitation spectrum by two-photon 
Bragg spectroscop y ^ 17 * 2 ^ employing amplitude modulation of the lattice po- 
tential. For the one-dimensional situation we modulate the weak axial po- 
tential, for the two-dimensional situation we modulate one of the two weak 
standing wave laser fields. The modulated lattice potential takes the form 
V y (y,t) = V y fi(l + A mo dSui{2TTv mo dt))sm 2 {ky). The modulation with am- 
plitude A mo d and frequency v m od generates two sidebands with frequencies 
iz^mod relative to the lattice laser frequency which define the energy hu mo d 
of the excitation. Due to the Bragg condition atoms scattering two photons 
receive a momentum transfer of Ohk or ±2hk. 

In contrast to applying a potential gradient across the lattice 1 ^, this 
method is not susceptible to effects like Bloch oscillations, dynamical in- 
stabilities, and Zener tunnelling which occur for low axial lattice depths. 
Furthermore, the excitation energy is precisely determined. 
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Fig. 2. a) Experimental sequence to perform Bragg spectroscopy on atoms 
in an optical lattice. The figure illustrates an experiment with a one- 
dimensional gas. For a two- or three-dimensional gas the curves are ac- 
cordingly different (see text), b) Scheme for resonant excitation of a Mott 
insulator by Bragg spectroscopy (not to scale). 




After the excitation, the experimental sequence continues by ramping 
down the lattice potentials linearly in 15 ms to V xy z = 4 -Br where the 
atoms are able to tunnel again in all three dimensions between the sites 
of the lattice. To allow for re-thermalization of the system, the atoms are 
kept at this lattice depth for 5 ms. Then all optical and magnetic potentials 
are suddenly switched off. The switch-off time for the optical lattice is 5 /xs 
and for the magnetic trap 300 fis. The resulting matter wave interference 
pattern is detected by absorption imaging after 25 ms of ballistic expansion. 
The width of the central momentum peak is taken as a measure of how 
much energy has been deposited in the sample by the excitation. If the 
energy increase is small, the peak is well fitted by a bimodal distribution. For 
resonant excitation there is only a single gaussian component, reflecting that 
the temperature of the atoms has significantly increased. To be independent 
of the shape of the peak we use the full width at half maximum (FWHM) 
as a measure of the introduced energy. Although this underestimates small 
energy increases, the important resonances and features of the spectra are 
well shown. 

We have chosen a duration t mo d = 30 ms and amplitude A mo d = 0.2 
of the modulation such that the resulting excitation of the condensate does 
not exhibit saturation effects for all measurements presented here (see Fig. 

We have experimentally verified that all atoms remain in the lowest 
Bloch band by adiabatically switching off the lattice potential^ after the 
modulation. For a thermal cloud of atoms loaded into an optical lattice we 
do not observe an energy increase by the modulation. 
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Modulation time t mod [s] Modulation amplitude A 

Fig. 3. Linearity of the excitation scheme for a one-dimensional gas. Filled 
symbols correspond to the superfluid regime, with V y = AEr, open symbols 
to the Mott insulating regime with V y = 10 Eft. All data where taken close 
to the peak of the energy absorption at a modulation frequency of v m od — 
1890 Hz. a) Modulation with constant modulation index A moc [ = 0.2. b) 
Modulation with fixed duration of 0.03 s. 



4.2. Results 

The fundamental change in the excitation spectrum when the system 
undergoes the quantum phase transition from a superfluid to a Mott in- 
sulator can be seen in Figure |1J The broad continuum of the superfluid 
(U I J <C 5.8) contrasts with the discrete spectrum of the Mott insulator 
(U/J > 5.8). Figure displays the ID situation with V± = 30 E R and 
includes a series of spectra for different values of U/J. The transition in 
a two-dimensional gas (V z = 30 Er) is shown in b). Finally, Figure 
represents the three-dimensional case with V x = V y = V z . One surprising 
feature is that we can excite the superfluid with our scheme at large hv mo( i 
contrary to predictions for the weakly interacting superfluid in an optical 
lattice formed by a single standing wave^. In our experiment strong inter- 
actions lead to a significant quantum depletion. It amounts to ~ 50% for the 
ID configuration with U/J = 2.3 3 ^. This has been studied quantitatively 
in a measurement of the coherence properties of the system^. Therefore 
this parameter may not be regarded small as in standard Bogoliubov theory 
and higher order excitations should be taken into accounlP^. In combina- 
tion with the broken translational invar iance in the inhomogeneous trap, 
this could explain the non-vanishing excitation probability observed in the 
experiment at high energies® . 

In the superfluid regime we obtain spectra which differ significantly 
from previous results for three dimensions^, since the superfluid excitations 
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Fig. 4. The measured excitation spectra of the superfluid to Mott insulator 
transition, a) An array of ID gases (V x = V z = 30 Er). b) Superfluid 
to Mott insulator transition in 2D (V z = 30 Er). c) Superfluid to Mott 
insulator transition in 3D (V x = V y = V z ). The ratios U/J given in brackets 
are calculated numerically. The data for a) and c) are taken from reference^. 
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decrease at higher energies. This decrease is rather slow for the ID gas 
but becomes more pronounced for the 2D and 3D gases. Our excitation 
scheme does not induce dephasing that occurs when the strongly interacting 
condensate is accelerated near the edge of the Brillouin zone^. This might 
cause the broadening and the background in the tilted lattice experiments^ 
at high energies. The width of the superfluid spectra for the ID gas is 
on the same order as twice the width of the lowest band for Bogoliubov 
excitations^. 

We observe the appearance of the discrete structure, which is character- 
istic for the Mott insulating phase, between U/J = 4 and U/J = 8. Above 
U / J = 20 there is no more background due to the vanishing superfluid com- 
ponent. However, the finite size and the inhomogeneity of the trap prohibit 
a sharp transition, so that the fraction of Mott insulating atoms increases 
gradually with U/J. In the Mott insulating phase we find the first resonant 
peak for all data sets close to the calculated value of U. A second peak ap- 
pears at (1.91 ± 0.04) times the energy of the first resonance. We attribute 
this resonance to defects where lattice sites with n = 1 atom next to sites 
with n = 2 atoms are being excited. For the ID and 2D system (Fig. 0^, 
and b) a much weaker resonance appears at (2.60 ±0.05) times the energy of 
the first resonance which could indicate higher order processes of two atoms 
tunnelling simultaneously. 



5. CONCLUSION 

We have studied the superfluid to Mott insulator transition in one, 
two, and three spatial dimensions. The excitation spectra were obtained by 
employing Bragg spectroscopy. We find the spectra for the low-dimensional 
superfluid to differ significantly from the three-dimensional case. 
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